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Abstract—A numerical model on the flexibility method in the case of a multilayer beam finite element has been developed 
and the contributions to its recent developments being made at Mechanical laboratory, Department of physics, Faculty of 
Sciences Rabat (Morocco). The results of the experiments and those of numerical calculations were concordant in the 
case of quasi-static loading. These results were based on the approach "finite element" coupled with a non-linear model 
[23]. Firstly, we present here the results based approach "finite element" related to the analysis of a bending square plate 
under concentrated and uniform load, clamped or simply supported on the contour. On the other hand, we present some 
results which we evidence to the problem related to the shear locking. The numerical model is based on a three- 
dimensional model of the structure seen here as a set of finite elements for multilayered plates multi cellular matrix 
(concrete) and a set of finite element fibers for reinforcement. The results obtained confirm the ability of these tools to 
correctly represent the behavior of quasi-statics of such a complex system and presage the deepening of a digital 
tool developed. 
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I. INTRODUCTION 

The phenomenon related to blockade by shear (or appearance of a parasitic stiffness) is a numerical problem that 
drew attention of many researchers in the past twenty years and an abundance of solutions which has been discussed in [3, 9, 
10, 11, 12, 19, 20, 22]. One way to avoid the appearance of shear locking and thus make the solution independent of the 
slenderness ratio (the ratio of length L I thickness h) is to calculate the terms of the stiffness matrix by integrating accurately 
the relative terms bending and sub-integrating the terms relating to shear [4,5,6,8,13,14,15,16,17,21 ,22]. To improve this 
phenomenon related to the numerical computation and propose a more efficient solution, we developed a model based on the 
flexibility method [23]. The model is formulated on the basis of the forces method by an exact interpolation stresses [18]. 
This makes it possible to calculate the flexibility matrix, which is the inverse of the stiffness matrix. The purpose of this 
study is the modeling of the structural response of the sails carriers subjected to seismic effects using a comprehensive three- 
dimensional numerical model using a nonlinear finite element approach coupled with a damage model developed for the 
behavior of concrete material. In this second paper, drawing on the results of the first article and those of [1,2 ,7], we present 
only some results related to the analysis of a homogeneous square plate in bending subjected to a concentrated and uniform 
load. 

II. MODELING 

Complementary to the trials and their interpretation, numerical modeling of this situation type has several 
advantages. In this case, it already developed an ambitious and effective model capable of taking into account the different 
aspects of this complicated problem, including the quasi-static and dynamic loading. Then after this satisfactory model, it 
has to constitute a way to complement the experimental measurements by providing new data. As such, it should contribute 
to a better understanding of the phenomena involved and to further provide a basis for dimensionality development methods. 

1. METHODOLOGY 

An immediate challenge before addressing the simulation of such problems is to choose the right methodology. 
The philosophy retained here is to realize the contribution of research in civil engineering to respond in a context of 
operational engineering. The choice was made on the use of finite element plate's multilayer multistage three nodes and two 
degrees of freedom per node. 

A realistic numerical prediction of the structural response of such a structure requires a rigorous three-dimensional 
geometric model of the system components. This model and its numerical analysis are implemented in the finite element 
code RE-FLEX. 

Then, the plate is meshed by including its geometry in a full mesh adapted to the different areas of the problem (it 
is discredited into layers and its thickness h in cells along x and y the surface) [Fig.l]. 



47 



Multicellular Multilayer Plate Model: Numerical Approach And Phenomenon . 




Figure 1 - Finite Element Model: Efforts resulting in a plate 

Where N xx ,N represent the normal forces and N the shear plane. M xx , M represent the bending moments and 
M torque. T x , T are the transverse shear stresses. 



2. CALCULATING THE ELEMENTARY FLEXIBILITY 

The exact interpolation functions are obtained by writing the various external forces of any point of the finite 
element, which here are the internal forces of the structure, according to the nodal reduced effort. Thus, we determine the 
matrices representing the exact interpolation functions of effort. The external forces of 'finite element' are supposedly similar 
with the same nature as the internal forces of the same element. 

One of the methods to calculate the external forces of "finite element" is the linearly interpolated from the equilibrium 
equations of the system. Notably in our study efforts are assumed constant at every point of "finite element" and moments 
vary linearly as a function of 

its variables (x and y in case of a plate). Thus, for a triangular plate finite element UK, we obtained the following 
relationships: 

The matrix that binds the membrane and bending efforts on any point with the reduced efforts is defined by: 

{*v} = {N,M} T = {N^N^N^M^M^Mj = [D^ij)]^} (i) 

The matrix that binds the shear efforts on any point with reduced efforts is defined by: 

{r} = {^r y } r = fc]{E;} (2) 
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[*-] = 



0 0 0 ^ 0 b 2 b 3 0 b 4 b 5 0 b t 
0 0 0 0 b 2 b x 0 b A b 3 0 b 6 b 5j 



(4) 



iJj iJj iJj jJj )Ji iJj 

5j = y t (X K -Xj)+ yj (Xj ~X K ) + y K (Xj -X,) is twice the area of the triangle UK 

{S; } T = { Nl , N 2 , N 3 , , M m , Mxyj , M XXj , M yyj , M xyj , M XXr , M„ k , M % } 

(5) 

Where j^j the vector of nodal efforts reduced, ^D cm j (£, T])~^ and [& c( ] are the matrices that represent accurate 

interpolation functions of the efforts membrane bending and shear respectively in the absence of apportionment. The 
stiffness matrix is simply the inverse of the flexibility matrix. 

|Z cm y | and | T x , T y | are respectively the vector normal forces, effort membrane, bending moments, twisting moment and 
shear forces applied to the cell. 

The direct connection of the finite element provides the stiffness matrix of elementary model in the local 
coordinate expressed by: 

[F>[*r[i?j>] 

Oil' =([F/>m/)] + [F/>«)]f (7) 

Where |^i^ a (c?n/ : ')^| and ^F^ (cis)~~^ are respectively the flexibilities of the matrices membrane combination 

bending and shearing of the plate. is the transition matrix to the system without rigid modes of deformation within five 

degrees of freedom, whose force field is represented by equation (8) and the corresponding displacements j^fj are defined 
(eqt.9): 

[F pla "] = [R] T {z;} (8) 
{q}=[R]{u} (9) 

With ^p plaq t ne external force exerted by a plate finite element nodal loads equivalent to the same element and jw e j the 
corresponding vector of nodal displacements and is given by equation (10): 

{u e } T = {u 0i , v 0i , w 0i , p xi , P yi , u Qj , v 0j , w 0j , P Xj , P yj u 0k , v 0k , w 0k , P Xk , P yK } (10) 

Remark: In the simple case of a beam with two nodes with three degrees of freedom [23] the force vector corresponds 
exactly to the demands of the nodal finite element beam. 

Flexibility matrices concerning the plates are given by: 

mcells 



[F£(cmfj] = S IJK £ JJ[z>^(^i7)] r [H cmf ^ K ,l K )Y [D cmf {$, V )\d$d V (ll) 

k=l A 
mcell s 

[F£(cisaM)] = S IJK £ \\[b ct ] T [HJ^tj,)]- [b^d&rj ( i 2) 
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The matrices cm f (<% K , JJ K ) J and \H ct (£ K , Tj K )] are matrices named flexibilities membrane bending and shear 
respectively: 
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The matrices ^H cm j(i^ K ,TJ K )~^ and [/f c (<f^, ^7^)] respectively represent the stiffness of membrane bending and 
shearing of the cell k of the plate, hi and Zi represent respectively the thickness and position Z layer i of the cell, E i and V i 
being respectively the Young's modulus and Poisson's ratio of the corresponding layer, k is the shear correction factor. 
\d C mf 0=^'%:)} = l^tt' e yy-> Y X y-> ^xr' ^yy ' /£iy } ' s tne vector °f plane deformation, and membrane of curvature 
experienced by a cell, and \d ct } = j^, / y | is the vector of deformations of the distortion in the planes (x, z) and (y, z). 

3. PRESENTATION OF AN ELEMENT DKT (Discrete Kirchhoff Triangle) 

The DKT element defined in [1] is a finite element with three nodes and three degrees of freedom per node. It is 
considered in this article, as a finite element with three nodes and five degrees of freedom per node. 

The rotations p x , /3 are interpolated in a parabolic manner and the transverse displacements U 0 , V Q , W Q are interpolated 
in a linear manner [1,2]: 

Px = j_>,A, + £ P«a k ,P y =t, N A + £ ■ P ^ = P ^ and p » = P A 

k=n+\ (=1 k=n+l 

n n n 
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Where C, , 5, are the direction cosines, k is the middle of respective sides of the triangle, and are given by the side 

if. C k = (Xj -x,)IL k , S k = {y j - y t )IL k and L k = J( Xj -xf +(y . -y t f 

Where n is the number of nodes of the finite element, in the case of a triangular element n = 3 and functions N i and P k 
are given by [1,2]: 

N l= A = l-^-r/, N 2 =%, N 3 =?j, P 4 = 4%A, P 5 =4^r/ and P 6 =4rjA 

The expression of CC k according to the nodal variables of nodes i and j is [1]: 
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a k =— (w, " Wj ) - ~ (Q A, + S k j3 Xi + C k j3 X/ + S k j3 yj ) 



(13) 
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III. ANALYSIS OF A UNIFORM PLATE WITH DKT AND FLEXIBILITY (FLX) 

At first glance, the figure 2 represents the results obtained with FLX as we analyze a homogeneous square plate 
subjected to uniform load simply supported or built on the contour, for different slenderness L/h (5 to 1000). The plate is 
meshed with 128 (N = 8) rectangular isosceles elements (§ 2.2). The results are virtually identical with those obtained with 
DST and Q4y [1] for the recessed plate (Figure 2. a). For the simply supported plate there appeared an error of 

about -0.5% (Figure.2.b). 
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b) Simply supported plate 
Figure 2 -homogeneous square plate with load uniform. Bending in the centre based on L/h 

(D = Eh 3 /12(l-v 2 ) ,v = 0.3, k =5/6) 
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In a second step, we describe in Figures 3 to 7 some results [1, 2] and on the analysis of a homogeneous square plate 
subjected to a concentrated load at the center or simply supported and built on the contour. A quarter of the plate is meshed 
with 2, 8, 32, 128 (N = 1, 2, 4, 8) rectangular isosceles DKT elements (§ 2.3) and FLX (§ 2.2). These elements have five 
degrees of freedom per node and are of Kirchhoff (no transverse shear energy, the results are independent of L/h) for DKT 

elements and flexible elements taking into account of transverse shear (FLX) for LI h = 24 (figures 6 and 7). Figures 4 
and 5 we presents the results obtained with FLX as we analyze a homogeneous square plate subjected to concentrated load 
simply supported or built on the contour, for different slenderness L/h (5 to 1000) for both types of mesh (there is a thin 
outlook of influence mesh ). The plate is meshed with 2, 8, 32 and 128 (N = 1,2, 4,8) rectangular isosceles FLX elements 
(§ 2.2). Convergence can be seen for N = 8, that is to say, for a mesh of 128 elements, we observe a occurrence of an error, 

for the clamped plate, in the order of —0.5% mesh A (Figure 4.a) and —0.56% mesh B (Figure 5. a) and simply 

supported plate —0.3% mesh A (Figure 4.b) and —0.31% mesh B (Figure 5.b). In Figure 6, we provide the percentage 
error of the deflection at the center depending upon 'N' number of divisions per half side. There is a monotonic convergence 

1 



with FLX (FLX model is a consistent shift, the total potential energy E p > E p 



and as E p — 



W. 



.P , we observe 



thatlw ) <(w ) ). It is observed that DKT is a model that over-estimates W . However, the monotonic 

V L ' EF V ' exact ' 

convergence of DKT can't be demonstrated. There is also a strong influence on the orientation of the mesh with triangular 
elements of the type DKT and FLX. The convergence of the moment | M x | in the middle of the recessed side and of the 

reaction concentrated in the corner (= YlM ) in case of simply supported plate are presented in Figure 7 for both types 

of meshes and for DKT and FLX, (Calculations of efforts have been made directly to the nodes peaks followed by an 
average if the node is shared by two elements). There is a fairly rapid convergence, an influence of models and an orientation 
of the mesh. 
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• Symmetry conditions: 

J3 X = 0 on CA ; on CD /3 y =Q 

• Boundary conditions: 

- Recess : w = P x = P y = 0 on ABD 

- Support simple: w = /3 X =0 on AB, 
w = p y = 0 on BD 



• Meshes considered: N=l, 2, 4,8 



Case N = 2 




mesh A mesh B 



Kirchhoff solution for a concentrated load P: 

(p = Eh 3 1 12(1 -v 2 );v = 0.3) 

Recess: W=5.6xl0" 3 PL 2 /D and \M | = 0.1257P 
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Simple Support: W c = 1 1.6 X 10 J PL I D and R = |2M vy | = 0.1219P 



Figure 3-square plate under concentrated load. Data 
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Figure 4-homogeneous square plate with concentrated load. Arrow in the center in terms of L/h (mesh A) 
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Figure 5 -homogeneous square plate with concentrated load. Arrow in the center in terms of L / h(mesh B) 



Where Wt the numerical value calculated for the different divisions (N = 1, 2, 4.8) 
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Figure 5-square plates with concentrated load at center built and simply supported. Error for DKT and FLX 
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IV. CONCLUSION 

The flexibility method developed with a linear interpolation (interpolation functions of the first order) and of way 
independently of the transverse displacements and rotations, solves the problem related to the phenomenon by 
shear locking. In the case of multicellular multilayer finite element, we observe that the method of flexibility, 
which is a model monotone convergence, converges quickly enough for a plate structure. In this paper we have 
presented the results for the analysis of a square plate in bending under load concentrated at the center, simply 
supported on the contour or clamped while highlighting the influence of the mesh on different slenderness L / h 

(Figures 4 and 5: arrow report W c I W k ). We also presented results on an analysis of a square plate subjected to a 

uniform load, clamped or simply supported on the contour (Figure 2). The percentage error appeared in Figures 4, 
5, 6 and 7 and that can be translated by the phenomenon of blocking is reduced (becomes negligible) by increasing 
the number of elements this allows us to confirm the reliability of the method on solving the problem of shear 
locking. In the following work (in a future article) we present the results at predictive calculation of the 
performance of bearing subject to the sails seismic behavior by numerical simulation coupled with a damage 
model by comparison with experimental results and by adopting a damage model for multicellular multilayer 
finite element . 
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